arXiv:1508.02866vl [astro-ph.HE] 12 Aug 2015 


High energy cosmic ray self-confinement close to extragalactic sonrces 

Pasquale Blasi^’^, Elena Amato ^ and Marta D’Angelo^ 

^INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi, 5 50125 Firenze, Italy and 
^Gran Sasso Science Institute (INFN), Viale F. Crispi 6, 60100 L’Aquila, Italy 

(Dated: August 13, 2015) 

The ultra-high energy cosmic rays observed at the Earth are most likely accelerated in extra- 
galactic sources. For the typical luminosities invoked for such sources, the electric current associated 
to the flux of cosmic rays that leave them is large. The associated plasma instabilities create magnetic 
fluctuations that can efficiently scatter particles. We argue that this phenomenon forces cosmic rays 
to be self-confined in the source proximity for energies E < Ecut, where Ecut « GeV for 

low background magnetic fields {Bo -C nG). For larger values of Bo, cosmic rays are confined close 
to their sources for energies E < Ecut ~ 2 x 10® GeV, where B-io is the field in units 

of 0.1 nG, Aio is its coherence lengths in units of 10 Mpc and L 44 is the source luminosity in units 
of 10"‘"‘ erg/s. 

PACS numbers: 98.70.Sa 


Introduction - The sources of ultra high energy cos- 
mic rays (UHECRs) are yet to be found. Actually at 
present there is no clear consensus even on the transition 
energy at which cosmic rays (CRs) start being mainly of 
extragalactic origin [l|. Some lower limits on the aver¬ 
age luminosity of UHECR sources have been presented 
in the literature @], hinting at the need for very high 
luminosities. 

In the following we will not be concerned about the 
nature of the source. We will only assume that UHECRs 
are accelerated in some unspecified extragalactic sources 
and focus on their propagation while leaving the parent 
galaxy. We will refer to the latter as the “source” for 
the sake of brevity, but the actual accelerator can be 
anything within the parent galaxy. 

Let us assume, for simplicity, that the CR sources ac¬ 
celerate a spectrum q{E) oc E~'^ up to some maximum 
energy Umax- The differential number density of CRs 
streaming out of such sources will read: 

« 1.7 X 10-1^7.44 cm-^GeV-i , (1) 

where we have taken A = ln(£’max/Einin) ~ 25 and 
Lcr = 10^"‘L44 erg/s, energies are in GeV and distances 
in Mpc. We assume that the source is in a region of 
the intergalactic medium (IGM), where the density of 
baryonic gas is ni, = ilbPcr/mp ~ 2.5 x 10“^ cm~^ 
(where pcr is the critical mass density of the universe 
and = Pb!Per — 0.022 is the ratio of the baryon den¬ 
sity pb with respect to the critical density) and we as¬ 
sume that there is a cosmological magnetic field with a 
strength Bq = G and a correlation scale Xb ~ 

10 Mpc Aio (where Aio is in units of 10 Mpc). Hence on 
scales smaller than Aio, the field can be considered as ori¬ 
ented along a given z direction. In such a situation, the 
Alfven speed is va = Bo/y/^'^^bPer ~ 44 i?_i 3 cm s“^. 


The positively charged CRs leaving their sources form 
an electric current that will be compensated by motions 
in the background plasma (return current) so as to en¬ 
sure local charge and current neutrality. This situation 
is known to give rise to a non-resonant plasma instability 
that is potentially very important for CR transport. This 
instability was first proposed in Ref. d, in connection 
with CR acceleration in supernova remnants, in which 
context it may provide a mechanism for strong magnetic 
field amplification in the shock proximity, a necessary 
condition to accelerate particles to ~ PeV maximum en¬ 
ergies. 

The main purpose of this letter is to evaluate under 
which conditions this phenomenon gives rise to an in¬ 
stability and what the consequences are in terms of CR 
propagation and intergalactic magnetic field generation. 


The calculations - The current associated with CRs 
streaming away from their sources in the IGM is easily 
written as a function of the minimum energy E of parti¬ 
cles in the current as 


JcK = encR(> E)c 


cLcrE ^ 

47rAr^ 47rr^ ’ 


( 2 ) 


where in the last equality we made use of Eq. [TJ 

This expression is strictly valid only if the background 
field is zero, or if the Larmor radius of the particles is 
CL 3> Xb, but we shall see that the above estimate for 
the current density turns out to hold also in the diffusive 
regime. 

A current propagating in a plasma can give rise 
to instabilities of different types. Granted that the 
current carrying particles are well magnetised {va > 
Jcr/( en-b), which from Eq. [5] is seen to imply Bq > 
2 X 10 T44Ag,gyr]y[p^) G, the fastest growing insta¬ 
bility arises when the condition 


JcrE > 


ccBq 

47r 


( 3 ) 
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is satisfied. This condition, which is the standard one for 
the development of non-resonant modes of the streaming 
instability Q, is equivalent to the requirement that the 
energy density of CRs be locally larger than the energy 
density in the form of pre-existing magnetic field, BqIAtt. 
For q{E) oc E~^, this requirement becomes independent 
of E and, using Eq. [2l can be simply formulated as: 

l1I2 

r < rinst = 3.7 x 10^—^ Mpc. (4) 

-D-13 

When Eq. S] is satisfied the fastest growing modes in the 
amplified field dB have a wavenumber fcmax that reflects 
the equilibrium between magnetic tension and Jcr x 6B 
force on the plasma, namely femax^o = ^JcK, and their 
growth rate is: 


Tma 


— ku 


kVa = 


I 47r Jcr 

nbirip c 


( 5 ) 


independent of the initial value of the local magnetic held 
Bq. The scale of the fastest growing modes is much 
smaller than the Larmor radius of the particles domi¬ 
nating the current (this is entailed in the condition for 
growth, Eq. El), therefore they have no direct influence on 
particle scattering. This conclusion is however changed 
by the non-linear evolution of the modes. As long as 
the instability develops on small scales, it cannot affect 
the current, hence one could treat the two as evolving 
separately. A fluid element will be subject to a force 
that is basically ~ JcrSB/c: its equation of motion is 

p{dv/dt) ~ \ Jcr5B, with 5B{t) = (5Bo exp (ymaxO- As 

an estimate, one can write the velocity of the fluid el¬ 
ement as u ~ i,^B{t)JcR)/icp"fma.x), which upon inte¬ 
gration leads to an estimate of the mean fluid displace¬ 
ment as Ax ~ JcR)(cp7n^ax)- We can then esti¬ 

mate the saturation of the instability by requiring that 
the displacement equals the Larmor radius of particles in 
the current as calculated in the amplified magnetic field, 
E/eSB: when this condition is fulfilled, scattering be¬ 
comes efficient and the current is destroyed. This simple 
criterion returns the condition: 


SB^ 

47r 


JcrE 


ce 


= ncR{> E)E. 


( 6 ) 


Since ncR(> E) oc E~^ in the case considered here, the 
saturation values of the magnetic field is independent of 
the energy of particles in the current driving the insta¬ 
bility. A somewhat lower value of the saturation was in¬ 
ferred in [^ , as due to the non-linear increase of the wave¬ 
length of the fastest growing modes. Following such a 
prescription our saturation magnetic field would be ^ 10 
times smaller. Eq. [6] expresses the condition of equipar- 
tition between the CR energy density and the amplified 
magnetic pressure, a condition that is often assumed in 
the literature without justification, and that here arises 
as a result of the physics of magnetic amplihcation itself. 


The field strength, as a function of the distance r will 
read 

5B{r) = 3.7 X Gauss . (7) 

This rather strong magnetic field will develop over dis¬ 
tances r from the source that satisfy Eq. 2] and under 
the additional condition that the growth is fast enough 
so as to reach saturation in a fraction of the age of the 
universe, to (numerical simulations of the instability Q 
show that saturation occurs when ^maxT ~ 5 — 10). This 
latter condition reads 7maxto ^ 5 and translates into: 

r < rgrowth = 1.2 X lO^Lli^E-%^ Mpc . (8) 

If the conditions expressed by Eqs. |4] and |8] are fulfilled, 
then the magnetic field can be estimated as in Eq. [7]and 
since SB ^ Bq and there is roughly equal power at all 
scales (equivalent to say that SB in Eq. [Slis independent 
of energy A), it is reasonable to assume that particle 
propagation can be described as diffusive, with a diffu¬ 
sion coefficient corresponding to Rohm diffusion in the 
magnetic field SB. This assumption is based on two dif¬ 
ferent considerations: Rohm diffusion regime is generally 
found in the quasi-linear theory of wave particle interac¬ 
tions when SB/Bq ^ 1 and SB^ is roughly independent 
on scale (the dependence is only logarithmic: see e.g. @); 
in the particular case of turbulence generated by Rell’s 
instability, with SB/Bq ^ 1 and scale invariant power 
spectrum, there is additional evidence that transport is 
governed by Rohm diffusion (see Q and 0 for exten¬ 
sive discussion). We can then write the particle diffusion 
coefficient as: 

D{E, r) = 9 X IO^^EgsV ^pc s'^ (9) 

The initial assumption of ballistic propagation of CRs 
escaping a source leads to conclude that particles would 
produce enough turbulence to make their motion dif¬ 
fusive. The diffusion time over a distance r from the 
source can be estimated as Td{E,r) = r'^/D{E,r) Ri 
3.3 X lO^^rMpc E/Zly from which follows that par¬ 

ticles can be confined within a distance r from the source 
for a time exceeding the age of the Universe, if their en¬ 
ergy satisfies the condition: 

E < Aconf = 2.6 X 10^ TMPC GeV . (10) 


One might argue that this conclusion contradicts the 
assumptions of our problem: for instance the density of 
particles in the diffusive regime is no longer as given in 
Eq.m This is certainly true, but the current that is re¬ 
sponsible for the excitation of the magnetic fluctuations 
remains the same, as can easily be demonstrated: for 
particles with energy E > Econf in Eq. [TOl and assuming 
that energy losses are negligible, quasi-stationary diffu¬ 
sion can be described by the equation 


J .2 Qj. 


r'^D{E,r)^ 


A'Kr'^ 


S(r), 


( 11 ) 
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where q{E) is the injection rate of particles with energy 
E at r = 0. Here the advection term has been neglected 
since there is no bulk motion of the background plasma 
and the Alfven speed is very small. This equation is 
easily solved to provide the density of CRs at distance r 
from the source: 


n{E, r) 


q 

SnrD^E, r) 


( 12 ) 


Clearly, by definition of diffusive regime, the density of 
particles returned by Eq. [T2]is larger than the density in 
the ballistic regime, Eq. [TJ However, the current in the 
diffusive regime is 


rdiff 


eED{E, r) 


dn 

dr 


47r7-2 


(13) 


which is exactly the same current that we used in the 
case of ballistic CR propagation (see Eq. [2|). This is a 
very important and general result: the magnetic field in 
Eq. [7] is reached outside a CR source independent of the 
mode of propagation of CRs, since it is only determined 
by the current and not by the absolute value of the CR 
density. Clearly the particles that are confined within a 
distance r around the source do not contribute to the CR 
current at larger distances. 


Results and implications - The confinement energy in 
Eq. [TU] is somewhat ambiguous since it depends on the 
distance r. What is the highest energy at which CRs 
escaping a source of given luminosity are confined to the 
source vicinity? In order to answer this question we need 
to take into account all the three conditions that need 
to be imposed to guarantee confinement, namely Eq. |4] 
(existence of fastly growing modes), Eq. [5] (growth rate 
faster than the expansion of the universe) and Eq. [10] 
(confinement). The first condition yields a limit on the 
distance from the source that is easy to satisfy, unless the 
strength of the background magnetic field is increased by 
several orders of magnitude, in which case however other 
complications arise (see discussion below). 

The other two conditions lead to the constrain 

Ecut ~ 10^ GeV L%^. (14) 

These particles are confined within a distance from the 
source: 


Mpc around the sources of UHECRs. As a consequence, 
the spectrum of CRs leaving these sources and eventu¬ 
ally reaching the Earth must have a low energy cutoff 
at an energy Ecut- This kind of cutoff has been often 
postulated in the literature in order to avoid some phe¬ 
nomenological complications that affect models for the 
origin of UHECRs. For instance, a low energy cutoff is 
required in the dip model [1, to describe appropriately 
the transition from Galactic to extragalactic CRs. This 
feature is usually justified by invoking some sort of mag¬ 
netic horizon in the case that propagation of UHECRs is 
diffusive in the lower energy part of the spectrum A 
similar low energy suppression of the CR flux is required 
by models with a mixed composition llj. In the calcu¬ 
lations illustrated above, the presence of nuclei is readily 
accounted for, provided the current is still produced by 
protons (assumed to be the most abundant specie). In 
this case, the value of Ecut is simply shifted to Z times 
higher energy for a nucleus of charge Z. 

As discussed above, the condition that guarantees the 
existence of non-resonant modes (Eq.|4|) is easily satisfied, 
unless the background magnetic field reaches Bq Ri 9.6 x 
10-10 . However, when this happens the calculations 

above break down for another reason: CRs can freely 
stream from the source only if their Larmor radius in 
the pre-existing magnetic field is much larger than the 
assumed coherence scale of the field, namely 


E > lOOCeV H-iaAio. 


(16) 


When Eq. |T 6 | is not fulfilled, namely when the back¬ 
ground field is relatively strong, then the propagation 
of CRs from the source becomes intrinsically one dimen¬ 
sional, which implies that the density of particles can be 
written as 


ncR{E,r) 


Q{E)t 

Trrj^vt 


2Q{E) 

■nRc{EYc 


(17) 


where we used the fact that the mean velocity of the 
particles is r; = c /2 for a distribution of particles that is 
isotropic on a half plane and we assumed that particles 
spread in the direction perpendicular to the background 
field by a distance equal to Rc{E) = max(ri(U), i?s) 
with Rg the source size and ri^[E) the Larmor radius 
of particles of given energy E. For a given source size 
tl > Rs as soon as Eq^y ^ 9 x 10®i?_io(i?s/100fcpc). At 
energies larger than this: 


Tconf ~ 3.8 Mpc l]^^. (15) 

Within such a distance the magnetic field is as given by 
Eq. 13 and larger than 5B{rcon{) ~ 9.6 x G. 

It is noteworthy that both the size of the confinement 
region and the magnetic field depend weakly upon the 
CR luminosity of the source, respectively as and 
L 4 / . Hence we can conclude that magnetic fields at the 
level of 0.1 — InC must be present in regions of a few 


ncR(A, r) ^ A7UiEQ%B\^cur^GeY-\ (18) 

If particles with energy > E are able to reach a given 
location, the current at that location is 

JcK - e.^EncK{E,r) = , (19) 

2 -nri 

which results in non-resonant growth of the field for 
E < Uinst = 3.5 X lO'^CeV 


(20) 
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independent of i?o, and in a saturation magnetic field 

SB « 0.7 B_io Gauss. (21) 

This value of the magnetic field is apparently very large 
and reflects the very large density of particles at low 
energies in the proximity of the source, as due to the 
reduced dimensionality of the problem. However one 
should notice that the value is normalized to the den¬ 
sity of GeV particles, which only live in the immediate 
vicinity of the source and generate small scale fields to 
which high energy particles are almost insensitive. At 
Mpc scales, where only high energy particles can reach, 
the field strength is much lower as we discuss below. 

Assuming again that the diffusion coefficient is Bohm- 
like, one can write: 

DiE, r) = 4.8 X 10^® E^^y cmVs, (22) 

which leads to an estimate of the diffusion time: rdifi = 
6.2 X L 44 ^ H_io r^pj, yr. Requiring that par¬ 

ticles reach the location at distance r in a time shorter 
than the age of the universe, we then obtain: 

rconf ~ 0.5 ^ lO^Gev) ^^ 44 ^ Mpc. (23) 

Following the usual procedure, one can calculate the 
growth rate of the fastest modes: 

^ s-\ 

(24) 

and impose the condition that 7 maxfo > 5, which reads: 

E < Agrowth ~ 5.3 X 10“ GeVL^^^H!./^ (25) 

The intersection of all the conditions listed above leads 
to conclude that particles with energies 

E < Flcut = 2.2 X 10®GeV B^J^^ Am (26) 

will be confined within a radius 

T-conf ~ 10 Mpc Aio . (27) 

The amplified magnetic field at such distance is 

5BK.ix 10-^G A^o^ (28) 

We emphasise again that the results illustrated both in 
the case of 3-d (lower Bq) or 1-d propagation (higher 
Bq) are only sensitive to the GR current, and hence in¬ 
sensitive to whether particle propagation is ballistic or 
diffusive. 

Discussion and Conclusions - It is often the case that 
CRs affect the environment in which they propagate, 
through the emission and absorption of waves that couple 


Tmax — 


ITk ecncR( 


them with the background plasma. The phenomenon of 
self-generation of waves is especially important close to 
shock fronts, where this process heavily affects the max¬ 
imum energy that can be reached [^, Il3 a nd is accom¬ 
panied by observational consequences [l3L Q, such as 
spatially thin rims of enhanced X-ray synchrotron emis¬ 
sion (see 15[ for a review). Here we investigated these 


processes when GRs start their journey from extragalac- 
tic sources: the escaping CRs form an electric current to 
which the background plasma, at density D^pcr, reacts, 
by generating a return current that in turn leads to the 
development of small scale instabilities. The growth of 
such instabilities leads to large turbulent magnetic fields 
and to enhanced particles’ scattering. 


The details of this process depend on the strength of 
the pre-existing magnetic field Bq'. if it is very weak (say 
< X0“40 then, in the absence of non-linear phenom¬ 
ena, CRs will try to propagate in approximately straight 
lines. The resulting electric current leads to the devel¬ 
opment of a Bell-like instability, that modifies the prop¬ 
agation of particles to be diffusive: particles with en¬ 
ergy < GeV are confined inside a distance of 

« 3.8Ly^^ Mpc from the source for times exceeding the 
age of the Universe, thereby introducing a low-energy 
cutoff at such energy in the spectrum of CRs reaching 
the Earth. Since the confinement distance is weakly de¬ 
pendent on the source luminosity, we conclude that a 
region with ~ nG fields should be present around any 
sufficiently powerful CR source. If larger background 
magnetic fields are present around the source, the gy¬ 
ration radius of the particles can be smaller than the 
coherence scale of the field, and in this case CR propa¬ 
gation develops in basically one spatial dimension. For a 
coherence scale of 10 Mpc, CRs are confined in the source 
proximity for energies E <2 x Aio GeV. 


It is currently not known whether the confinement phe¬ 
nomenon occurs in one or the other regime since only 
limits exist on cosmological magnetic fields: upper lim¬ 
its can be obtained from Faraday rotation measures 161 
but these limits are rather weak (< nG) and model de¬ 
pendent. A lower limit can be found from gamma ray 
observations of distant TeV sources [l3, [13] and these 
limits are typically at the level Bq > 10“^^ G. Numerical 
simulations of large scale structure formation in the pres¬ 
ence of background magnetic fields typically find ~ 10 “^^ 
Gauss magnetic fields in voids [l3| (although see [l^l for 
different conclusions). 


The physical prescription adopted here leads to esti¬ 
mating the strength of the self-generated magnetic field 
5B at the level of equipartition with the energy in the 
form of escaping cosmic rays, independent of the value of 
the pre-existing field, Bq. A weak dependence on Bq was 
instead found for the saturation level in [ 3 | , which in our 
case would lead to 6B about ~ 10 times smaller for small 
values of Bq, thereby reducing the energy below which 
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CRs are confined in the source proximity. Understand¬ 
ing the dynamics of the magnetic field amplification and 
saturation is clearly very important. One could test the 
amplification mechanism in the case of supernova rem¬ 
nant shocks: in this case the saturation criterion used 
here translates to SB Ri yjA ttwc'rvs/ c, with vs the ve¬ 
locity of the supernova blast wave and wcr is the energy 
density in accelerated particles. Applying this criterion, 
we obtain an estimate of the magnetic field which is in 
good agreement with that measured in young galactic 
SNRs [1^. On the other hand, due to the relatively small 
value of SB/Bq, the saturation provided by Q would re¬ 
turn a value oi SB only a factor ~ 2 smaller, too small 
a difference to discriminate between the two estimates. 
The testing is then left to numerical experiments study¬ 
ing the propagation of a current of energetic particles in 
a low density, low magnetic field plasma: hybrid simula¬ 
tions with this aim are currently ongoing. 

The phenomenon of CR confinement illustrated here 
has profound implications for the description of the tran¬ 
sition region between Galactic and extra-galactic CRs 
BiEl- It is rather tantalising that the cutoff obtained 
here as due to self-trapping is in the same range of values 
that have previously been invoked in the literature based 
upon phenomenological considerations. 
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